LAYER-BY-LAYER PATTERN PROPAGTION 
AND PULSED LASER DEPOSITION 

F. WESTERHOFF, L. BRENDEL, D.E. WOLF 
Gerhard- Mercator-Universitat Duisburg, 47048 Duisburg 

In this article kinetic Monte Carlo simulations for molecular beam epitaxy (MBE) 
and pulsed laser depositon (PLD) are compared. It will be shown that an optimal 
pattern conservation during MBE is achieved for a specific ratio of diffusion to 
deposition rate. Further on pulsed laser deposition is presented as an alternative 
way to control layer by layer growth. First results concerning the island density in 
the submonolaycr regime are shown. 



1 Introduction 

Multi-layered structures with sharp interfaces are of high interest, due to the 
new physical properties they exhibit resulting from the reduced dimensional- 
ity. In order to obtain such sharp interfaces the deposited material has to grow 
in layers instead of forming 3D islands. Two different methods of growing thin 
films will be compared here. 

Molecular beam epitaxy (MBE) allows to produce mono-crystalline layer 
structures with high sharpness at the interfaces. Growth takes place in an 
ultra high vacuum chamber in which the material is thermally evaporated 
and thus is deposited on the substrate with a continuous flux F. 

The fabrication of thin films with pulsed laser deposition (PLD) is a quite 
new technique which can improve layer by layer growth. 1 In addition PLD 
provides the possibility to deposit any desired material with the stoichiometry 
of the ablated material being conserved. In contrast to thermal evaporation 
in MBE many particles are ablated from the target simultaneously, resulting 
in a high adatom density directly after the pulse. This causes an increased 
islands density in the PLD as will be described. 

2 Molecular Beam Epitaxy 

For an initially flat surface layer by layer growth usually improves, the higher 
the temperature of the substrate is. 2 But a pattern (e.g. a structure produced 
by lithographic techniques) will in general not be conserved better during 
layer by layer growth, if the substrate temperature is raised. The reason is, 
that for high temperatures the adatoms do not nucleate on top of the islands 
any more because the diffusion length becomes too large. Instead they are 



incorporated at the edges of the islands and therefore flatten the surface which 
means that information about the pattern gets lost. As the diffusion length 
depends on the ratio of the diffusion rate to the deposition rate, 3 this value 
has to be optimized if a given pattern should be conserved during growth. 
One can also consider the opposite optimization problem: How to choose the 
growth parameters in order to heal any unwanted structures as efficiently as 
possible? As the two optimization problems are obviously two sides of the 
same medal, we consider only the first one in the following. 

2.1 Model 

In the simulation we use a simple cubic crystal with periodic boundary condi- 
tions and assume that growth takes place without overhangs and defects. The 
size of the substrate is (L ■ a) 2 , L 2 denoting the number of lattice sites and a 
the lattice constant. The atoms are deposited with a flux F and diffuse with 
the diffusion constant D until they are irreversibly bound to another adatom 
or an island edge. 

By pattern propagation we mean the deterministic height increase by one 
lattice constant during the deposition of one monolayer. The probability for 
pattern propagation 4 is then the fraction of sites which were conserved in 
every monolayer until time t: 



2.2 Results 

Figure 1 (a) shows p(t) for an initially flat surface. The exponential decay 
p(t) = cxp(—t/t c ) is what one would expect, if the propagation of the surface 
sites from one layer to the next would be temporally uncorrelated with con- 
stant transition probability p(l). Surprisingly, the exponential decay of p(t) is 
found although there is clear evidence that transition probabilities are tempo- 
rally correlated. 5 The life time t c (Fig. 1 b) of the initial surface configuration 
is determined from the slope of ln(p(t)). The data give a power law for t c 
extending over 6 decades: 



In the limit of D/F — > (random deposition) we get t c = l. 6 That the last 
data point in Fig.l (b) deviates from the power law is a finite size effect. 
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Figure 1. (a) Propagation for an initially flat surface, (b) Life time t c (D/F). 



In a study of the one dimensional case 4,6 a logarithmic dependence 
(t c ~ \og(D/F)) was found. Most recent calculations seem to show a cross 
over from power law to logarithmic dependence for larger D/F in two dimen- 
sions, as well. 7 

2.3 Structured surfaces 

As long as the pattern is large compared to the diffusion length for an initially 
flat surface, the propagation probability should not change since the adatoms 
will not notice the pattern during their life time. Therefore we also expect an 
exponential decay for small ratios of D / F. Here we consider larger values for 
D/F, though, where we observe a different behavior. We find a maximum in 
t c (D/F) which is shifted to higher D/F with growing pattern size r (Fig. 2 a). 
This can be explained by the existence of two different mechanisms responsible 
for the pattern's decay. 

For small D/F the life time t c is independent of the island size. This is 
the regime where the pattern vanishes through kinetic roughening. To the 
right of the maximum the diffusion length Id is large compared to the pattern 
size, hence there are almost no nucleations on top of the "artificial" islands. 
Instead, they take on the form of natural growing and coalescing islands. In 
this way the pattern is also lost. 

We find (see Fig. 2 b) that the maximum life time t Cima x for a given pattern 
length scale r is reached for 

(D/F) max oc r 2 - 4 ^' 1 andscaleshke t c , max oc r°- 50±0 - 05 . (3) 
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Figure 2. (a) Life time t c for patterned surfaces, (b) Scaled data. 



This is in agreement with (2): i c , max cx {D/Ff m % 2A = {D/F)%£. As (2) 
crosses over to a much weaker dependence for very large D/F, we expect also 
that (3) will be modified for larger r. 

3 Pulsed laser deposition 

In pulsed laser deposition the atoms are deposited with a pulse intensity / 
(number of atoms per unit area) and frequency u> forming an effective (aver- 
age) flux 

F = I-u. (4) 

If we reduce the intensity / to the limiting case of one atom per pulse 
and at the same time increase the frequency so that the effective flux stays 
constant, we will reach conditions equivalent to thermal deposition (MBE). a 
One should therefore expect that both methods show similar characteristics 
in the submonolayer regime. 

If the intensity is increased, one should expect a higher island density, 
because the adatom density and therefore the probability for nucleation is 
also larger. Hence we predict a critical pulse intensity I c which separates 
both regimes. 

In the simulation the kinetic energy of the arriving atoms is neglected. 
This should be a valid approximation for experiments, in which the laser 

"Nevertheless in the simulation there is still a small difference since in PLD atoms are 
deposited at fixed time intervalls, whereas in MBE fluctuations may occur. 
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intensity is just above the ablation threshold 1 or the deposition takes place, 
e.g. in an oxygen atmosphere which slows down the arriving particles. 

3. 1 Critical pulse intensity 

During thermal deposition the time evolution of the adatom density p is given 
by a gain term, the flux F and a loss term p/r, where r is the life time of 
adatoms on the surface: 

p = F- 9 -. (5) 

T 

In the quasistationary regime (after island nucleation and before island coa- 
lescence) the life time r is determined by the typical distance Id an adatom 
diffuses before reaching an island edge 8 : t = Id 2 /D. Then p can be neglected 
and one gets: 

p^F-r^F 1 ^. (6) 

With Id oc (7^) 7 and 7=1/6 for two dimensions and irrevsible aggregation 3 
this leads to 

pcx/zT 4 . (7) 

Equation (7) describes the adatom density for thermal deposition. In PLD 
with a pulse intensity / much smaller than this density, the adatom density will 
oscillate around this value with the frequency u). Then one expects the same 
scaling of the island density as for thermal deposition. By contrast a much 
larger pulse intensity would prevent the adatom density in PLD to become 
on average similar to the one in MBE. Hence the critical pulse intensity is 
expected to be of the order 

h oc l D - A cx - . (8) 

In the case of a large pulse intensity, the island density increases to such a 
degree that the life time of the adatoms is shorter than the time between two 
pulses. Hence the adatom density and therefore the island density becomes 
independent of the pulse frequency and only depends on /. 

3.2 Results 

Fig. 3 (a) depicts the island density as a function of the deposited monolayers. 
We expect three different regimes 9 : The regime for the nucleation of islands, 
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Figure 3. (a) Island distance l]j as a function of time t (in units of monolayers), (b) 
Id{D/F) for the limiting case of one atom per pulse. 



the saturation of the island density and finally the coalescence of islands with 
an increasing island distance. An increase of the pulse intensity leads to an 
overall decrease of the island distance. 



3.3 Id for the limit of small pulse intensities 

According to a classical result, confirmed by experiments 10 and simulations 9 ' 11 
the island distance scales like: 

Id oc i^-j (9) 

This behavior is confirmed in figure 3 (b) for the limit, that only one atom 
per pulse (/ • L 2 — > 1) is deposited. 



3-4 Scaling behavior of the island distance 

Figure 4 (a) depicts the island distance Id as a function of the pulse intensity /. 
For small values of / the diffusion length approaches a constant value, which 
is the island distance for the MBE case. For larger pulse intensities Id is 
described by a power law: 

l D ozr v . (10) 
It turns out that Id can be written in a scaling form: 
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Figure 4. (a) Island distance as a function of the pulse intensity for different ratios of D/F. 
(b) Scaled data. 



with 

f = const, for y« land 

1 ~ 2/ for y > 1 . y ' 

While the asymptotic power laws (9,10) were also obtained by Jensen and 
Niemeyer 12 , the scaling law (11) is new. In the previous work the pulses 
had a finite duration in contrast to the case presented here. This additional 
characteristic time spoiled the scaling. 

Since for large /, the w-dependence (via F (4)) has to vanish: 



D 
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^•[p] = const - (13) 



and therefore: 

h oc ^ j (14) 

For large pulse intensities we find in Fig. 4 (a) the slope v — 0.248, so that, 
together with 7 = 1/6, we get: 

7/1/ w 2/3, 

in accordance with the theoretical prediction (8). 
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